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ABSTRACT 

K-dwarfs have life-times older than the present age of the Galactic disc, and are 
thus ideal stars to investigate the disc's chemical evolution. We have developed several 
photometric metallicity indicators for K dwarfs, based an a sample of accurate spectro- 
scopic metallicities for 34 disc and halo G and K dwarfs. The photometric metallicities 
lead us to develop a metallicity index for K dwarfs based only on their position in the 
colour absolute-magnitude diagram. Metallicities have been determined for 431 single 
K dwarfs drawn from the Hipparcos catalog, selecting the stars by absolute magnitude 
and removing multiple systems. The sample is essentially a complete reckoning of the 
metal content in nearby K dwarfs. We use stellar isochrones to mark the stars by mass, 
and select a subset of 220 of the stars which is complete in a narrow mass interval. 
We fit the data with a model of the chemical evolution of the Solar cylinder. We find 
that only a modest cosmic scatter is required to fit our age metallicity relation. The 
model assumes two main infall episodes for the formation of the halo-thick disc and 
thin disc respectively. The new data confirms that the solar neighbourhood formed on 
a long timescale of order 7 Gyr. 

Key words: Stars - K-dwarfs, abundances; Photometry - Johnson-Cousin, Stromgren 
and Geneva systems 



1 INTRODUCTION 

A central issue in studies of the chemical evolution of the 
Galactic disc is to resolve the so called "G-dwarf problem" 
(van den Bergh, 1962; Schmidt, 1963; Pagel and Patchett, 
1975). The problem is that the observed stellar metallic- 
ity distribution shows far fewer metal deficient stars than 
the predictions of the simplest, closed box models of the 
Galactic disc's chemical evolution. Integrated light studies 
indicate that the G-dwarf problem is not restricted to our 
own disc (Worthey et al, 1996) but is also found in other 
galaxies. There are many ways that the evolutionary mod- 
els can be modified to bring them into consistency with the 
observations, such as pre-enrichment of the gas, a time de- 
pendent Initial Mass Function or gas infall (for a review see 
e.g. Pagel, 1997) 

G-dwarfs are sufficiently massive that some of them 
have begun to evolve away from the main sequence, and 
these evolutionary corrections must be taken into account 
when determining their space densities and metallicities. 
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While these problems are by no means intractable, it has 
been long recognized that K dwarfs would make for a cleaner 
sample of the local metal abundance distribution, because 
for these stars the evolutionary corrections are negligible. 

K dwarfs are of course intrinsically fainter, and it has 
not been until recently that accurate spectroscopic K dwarf 
abundance analyses have become available, with which 
to calibrate photometric abundance indicators (Flynn and 
Morell, 1997). As a result, it is now known that there is a K- 
dwarf problem which is very similar to the G-dwarf problem 
(Flynn and Morell, 1997). Studies of M dwarfs indicate that 
the problem is present in these stars too (Mould 1976). Such 
stars are still not really ideal for measuring the local metal- 
licity distribution, because metallicities for M dwarfs are dif- 
ficult to determine in the optical (but appears possible using 
infrared JHK photometry, (see e.g. Stauffer and Hartmann, 
1986; Leggett et al, 2000). The all sky surveys presently un- 
derway in the infrared (Denis, 2MASS) may make such stars 
viable metallicity tracers in the near future. 

The G-dwarf metallicity distribution is already an 
extensively studied subject (Pagel and Patchett, 1975; 
Sommer-Larsen, 1991; Wyse and Gilmore, 1995; Rocha- 
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Pinto and Maciel, 1996) and has been used by chemical evo- 
lution models to constrain the time scale of the formation of 
the Galactic disc at the solar neighbourhood. The G-dwarf 
metallicity distribution can be well fit by models in which 
gas has been settling onto the disc over a protracted period, 
of some billion years. For example, Chiappini et al (1997), 
on the basis of the fit of the G-dwarf metallicity distribu- 
tion, have shown that the disc in the solar neighbourhood 
was formed on a long time scale of 7 - 8 Gyr. This conclusion 
has been later stressed also by other authors (e.g. Portinari 
et al, 1998; Prantzos and Silk, 1998; Chang et al, 1999). 

The release of the data from the European Space 
Agency's Hipparcos (ESA, 1997) satellite offers a great op- 
portunity to determine the local metallicity distribution of 
the disc from a complete sample of K dwarfs. There are 
three clear improvements because of Hipparcos. Firstly, the 
parallaxes are so accurate that the K dwarfs can be selected 
by absolute magnitude rather than colour, which is a much 
better way of isolating stars in a particular mass range. Sec- 
ondly, the uniformity of the Hipparcos data allows us to 
construct samples with well understood completeness lim- 
its. Thirdly, a large fraction of the close binaries can be 
identified using Hipparcos and removed (since photometric 
abundance indicators are calibrated for single stars). This 
last effect turns out to be quite important. 

In this paper we develop a number of simple photo- 
metric metallicity indicators for K dwarfs, based on a spec- 
troscopically determined sample of metallicities by Flynn 
and Morell (1997). The metallicities for several hundred K 
dwarfs drawn from the Hipparcos catalog have been mea- 
sured via new observations described here. This sample has 
been used to show that there is a very simple relation be- 
tween the V band luminosity Mv of K dwarfs at a given 
B — V colour and metallicity (the work is fully described 
in a companion paper (Kotoneva, Flynn and Jimenez, 2002, 
paper II). This simple relation is used in this paper to deter- 
mine metallicities for 431 single K dwarfs in a near complete 
sample selected from the Hipparcos catalog. We examine a 
set of isochrones and find a simple relation between the po- 
sition of a star in the Hipparcos colour magnitude diagram 
and mass. We are thus able to select a sample of dwarf stars 
by the more physically relevant parameter of mass, rather 
than by colour or absolute magnitude, as is usually the case. 
The metallicity distribution obtained is then compared to 
the results of models of the chemical evolution of the Solar 
neighbourhood. 

The paper is organized as follows: In section 2 we de- 
scribe the selection of the K dwarf sample and in section 3 
the observations and reductions of the stars in broadband 
and intermediate band photometric systems. The details of 
the various metallicity calibrations are presented in section 
4. In section 5 we use isochrones to fit the masses of the stars, 
and select the stars in an appropriate mass range which is 
near complete; we discuss the kinematics of the stars and 
compute corrections for the metallicity distribution function 
from the solar volume to the solar cylinder. The metallicity 
distribution is compared briefly with other determinations 
in the literature. We compare the data to the predictions of 
a model of the chemical evolution of the Solar neighbour- 
hood due to Chiappini et al (2001) in section 6. We draw 
our conclusions and summarize in section 7. 



2 THE SAMPLE 
2.1 Sample selection 

Our sample of K dwarfs is drawn from the ESA Hippar- 
cos catalog (ESA, 1997). For the purposes of this paper we 
term K dwarfs to be stars of absolute magnitude My in 
the range 5.5 < Mv < 7.3. The choice of these absolute 
magnitude limits is motivated as follows. The upper mag- 
nitude limit at Mv = 5.5 is chosen to avoid the effects of 
stellar evolution. Examination of theoretical isochrones (see 
e.g. Jimenez, Flynn and Kotoneva, 1998) indicate that the 
effects of stellar evolution on luminosity at My — 5.5 during 
the disc lifetime amount to at most 0.1 magnitude, typically 
much smaller. The effects of stellar evolution in constructing 
the sample are thus small, and negligible compared to the 
main source of error (which is Poisson sampling statistics). 
The limit at Mv = 5.5 corresponds to a spectral type of 
about G8. The lower absolute magnitude limit at My — 7.3 
is the magnitude of the reddest K dwarfs for which our pho- 
tometric metallicity indicators can currently be calibrated 
via spectroscopic observations (Flynn and Morell, 1997). 
The limit at Mv = 7.3 corresponds a spectral type of about 
K5. 

Since we are interested in obtaining a complete sam- 
ple of K dwarfs in the Solar neighbourhood, with which 
to construct the metallicity distribution function (MDF), 
the stars were initially selected from the "survey" part of 
the Hipparcos catalog, which is complete to an apparent vi- 
sual magnitude given by V < 7.3 + l.lsin|fe|. Here b is the 
Galactic latitude (the apparent visual magnitude limit was 
made dependent on b in order to avoid observing excessive 
numbers of stars in the Galactic plane). Adopting this ap- 
parent magnitude limit and the absolute magnitude range 
5.5 < Mv < 7.3 resulted in a sample of 209 stars. Our in- 
tention had been to get a sample of order 500-750 stars. 
In order to increase our basic sample size, we increased the 
magnitude limit by 0.9 mag (i.e. by taking an apparent mag- 
nitude limit of V < 7.3 + 1.1 sin|6| + 0.9) we obtained 668 
stars (including the first 209 stars). This is the basic sample 
which we targeted for observations. 

We need to understand the completeness of this basic 
sample, since it is drawn from stars 0.9 magnitudes fainter 
than the magnitude limit of the complete part of the Hippar- 
cos catalog. Even 0.9 magnitudes beyond the completeness 
limit, the sample turns out to be still satisfactorily close 
to complete for our purposes. The completeness level is ap- 
proximately 94%. We determined this quantity by using the 
Galactic structure model of Holmberg (2000). The model is 
based on star count data and makes use of both the Hippar- 
cos and Tycho data (the latter is complete to much fainter 
limits than Hipparcos) to construct the luminosity function 
of the local disc. Using the model, we find that for an ap- 
parent magnitude limit of V < 7.3 + 1.1 sin|£>| + 0.9 and 
absolute magnitude limits of 5.5 < Mv < 7.3, we would 
expect that some 710 stars should be in our basic sample, 
whereas there are actually 668 such stars in the Hippar- 
cos catalog. We conclude that even 0.9 magnitudes beyond 
the "survey" completeness limit, Hipparcos is substantially 
(94%) complete in the absolute magnitude range of interest. 

Stars were observed by Hipparcos up to 4 magnitudes 
fainter than the faintest stars in the "survey". The 0.9 mag- 
nitude extension beyond the survey limit results in a sam- 
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Figure 1. The upper panel shows colour-magnitude diagram for 
431 stars in our sample for which B — V colour is available, while 
the lower panel shows the colour-magnitude diagram for the R — I 
colour. Note that the main sequence is clearly thicker as a function 
of B — V compared to R — I. This is a direct consequence of the 
metallicity of the stars, which affects B — V much more than it 
affects R - I (Flynn and Morell, 1997). 

pie which is still over 90% complete. Beyond this 0.9 mag 
limit, completeness drops rapidly, and the stars in Hippar- 
cos become increasingly dominated by objects included for a 
particular astrophysical interest, such as having a low metal- 
licity and/or high velocity. There will be a very small excess 
of such stars in our extended sample, but the high com- 
pleteness (94%) means that it is other uncertainties (Poisson 
statistics) which dominate the construction of the metallic- 
ity distribution. 

2.2 Removal of multiple stars 

The photometric metallicity indicators described in section 
4 were calibrated using spectroscopically determined metal- 
licities for single stars. We examined the effect that multiple 
stars would have on our metallicity indicators, by making 
Monte-Carlo simulations in which we combined the fluxes 
of a range of pairings of single K and M dwarfs selected 
at random from the Hipparcos delineated main sequence, 
and computing the effect on the photometrically determined 
metallicity. Metallicities for multiple stars were found to be 
as much as 0.4 dex lower than the true metallicity. Clean- 
ing the sample of multiple stars thus turned out to be very 
important. 

The Hipparcos catalog included a flag for "probable 
multiple stars", based on the "reliability of the double or 
multiple star solution". This flag was used to eliminate all 
definite, possible and suspected multiple systems. This re- 
duced the sample from 668 to 449 stars or about 2/3 of the 



initial sample. Despite this expedient, a small number of bi- 
naries seem to remain in the sample (of order 10%). This is 
evident from the positions of the suspected multiples in the 
colour-magnitude and two colour B — V versus R — I dia- 
grams. This issue is studied in detail in a companion paper 
(Kotoneva, Flynn and Jimenez, 2002, Paper If). About half 
of these extra suspected multiples could be removed in con- 
structing the final sample, so the final sample should have 
a minor contamination by multiples of less than 5%. Note 
that, by "single stars", we mean only stars which have no 
companion bright enough to significantly affect the metallic- 
ity measurement (i.e. within 5 magnitudes of the brightness 
of the primary). After removing these suspected multiple 
stars the final sample consists of 431 stars. 

The colour-magnitude diagrams for our final sample of 
single stars are shown in Fig 111 



3 OBSERVATIONS AND REDUCTIONS 

Observations were made at Siding Spring Observatory 
(SSO) in Australia, between 8 th and 22 nd of March 1999. 
We measured Johnson-Cousins and Stromgren colours for 
all available Southern hemisphere stars of the basic sam- 
ple using the SSO 24" telescope. We used the Motorized 
Filter Box with all eight filters U, B ,V, R, I and Stromgren 
v, b and y. UBVRI primary standards were selected from 
Landolt (1983a, 1983b, 1992) and among E-region standard 
stars (Graham, 1982; Menzies et al., 1989). Stromgren pri- 
mary standards were selected from the E-regions and also 
from Gr0nbech, Olsen, and Stromgren (1976) and Crawford 
and Barnes (1970). Hauck and Mermilliod (1998) stars were 
used as secondary standards for the Stromgren system and 
Flynn and Morell (1997) stars for both of the intermedi- 
ate band systems. We observed 218 stars and 55 of these 
more than once. The standard error in the R — I colour was 
ps 0.012 mag, quite accurate enough for our purposes. In the 
Stromgren bands the photometric error was slightly higher 
but still « 0.02 mag. 

Further data were obtained at La Palma using the 
Swedish 60 cm telescope in March 2000. An SBIG CCD was 
used to observe in the filters V, R and I. 90 stars were ob- 
served, most of them more than once. The scatter for the La 
Palma observations was bit higher than at SSO, the scatter 
in the R — I band being 0.021 for a single observation. Since 
most of the stars were observed twice the mean error came 
down to w 0.015 mag. 22 of the stars were observed both in 
Australia and La Palma; the mean difference in R — I was 
less than 0.01 mag and the scatter was 0.01 mag, so that no 
significant difference was found between observations taken 
at the different observatories. 

We also obtained R—I colours for 108 stars from Bessell 
(1990), who observed all the stars in the Gliese catalog. 
These are found to be in excellent agreement with our own 
SSO data with a mean difference in R — I of less than 0.01 
mag with a scatter of 0.013 mag. A comparison of the two 
sources of R — I colour is shown in Fig 0. 
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Figure 2. Comparison of our R — I colours obtained at SSO 
versus those from BesselFs (1990) values from his observations of 
stars in the Gliese catalog. The mean difference is less than 0.01 
mag and the scatter is 0.013 mag. 

4 DETERMINATION OF METALLICITIES 

We have obtained metallicities for all 431 stars of the 
dataset. In the end, a single method was used, although 
we discuss here four metallicity indicators for the dwarfs. 
Firstly, we used a relation based on Geneva and R — I pho- 
tometry due to Flynn and Morell (1997). Secondly, we have 
developed a very similar method to Flynn and Morell's based 
on Stromgren and R — I colours. Thirdly, we have found a 
simple relation between the broadband B — V and V — I 
colours and metallicity. Fourthly, we have developed in a 
companion paper (Kotoneva, Flynn and Jimenez, 2002) a 
metallicity indicator for K dwarfs based on their absolute 
magnitude in the V band relative to a fiducial solar metal- 
licity isochrone. The first three methods have a similar error 
0.2 dex in [Fe/H], while the last method appears to be bet- 
ter, with an estimated error of 0.1 dex. This last method 
depends on the availability of accurate parallaxes for the 
stars, and so is less generally applicable than those which 
rely on photmetric colours alone. The last method is the 
one we adopt here for measuring the K dwarf metallicities. 

4.1 Metallicities from Geneva photometry 

To obtain the metallicities for the K dwarfs for which a 
Geneva bi colour is available, we used existing relations de- 
scribed in Flynn and Morell (1997). The b\ colours were 
obtained from Rufener (1989) for 245 stars and the R — I 
colours come from our observations and/or from Bessell 
(1990). Both of the colours were available for 149 stars for 
which the metallicities, [Fe/H]^, were computed using the 
relation: 

[Fe/H] 6i = 8.248 x bi - 12.822 x (R - I) - 4.822. (1) 
These metallicities have a typical error of 0.2 dex. 
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Figure 3. Relation between the metallicities using Stromgren mi 
colour and spectroscopically determined metallicities. The scatter 
around the 1:1 line is f» 0.20 dex. 

4.2 Metallicities from Stromgren photometry 

We have developed a new metallicity indicator for the K 
dwarfs for which Stromgren colours are available. The cal- 
ibration was obtained using the 34 G and K dwarfs from 
Flynn and Morell (1997). For these stars accurate, spectro- 
scopically determined metallicities, [Fe/H] spec and effective 
temperatures have been determined with errors of 0.05 dex 
and « 100 K, respectively. 

We searched for a relation between R — I 
colour, Stromgren colours and spectroscopical abundance, 
[Fe/H] spcc . The best fitting Stromgren colour was found to 
be mi, with a small dependence on c\. The b — y colour was 
found to be poorly correlated with metallicity as one might 
expect. The relation we found is: 

[Fe/H] mi = 6.076 x mi - 3.811 x ci + 9.4303 x ci 2 (2) 
-12.822 x (R- I) +3.226. 

Fig ^ shows the relation between the spectroscopically 
determined metallicity and the Stromgren and R — I based 
metallicity. The scatter around the one-to-one relation is 
» 0.20 dex. 

For 142 stars both the Stromgren mi and ci colours and 
R — I colour were available and the metallicities [Fe/H] mi 
were calculated. The dependence on ci is quite weak, so that 
while for 54 stars no ci colour was available, for these stars 
we adopted a typical value of ci = 0.29. Adopting this mean 
value leads to a very small increase in the metallicity error 
of » 0.02 dex. 



4-2.1 A check using the Hyades 

A check of the Stromgren calibration was made using G 
and K dwarfs in the Hyades cluster (Reid, 1993; Flynn and 
Morell, 1997), using VRI photometry and Stromgren mi 
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Table 1. Broadband and Stromgren m\ data for G and K dwarfs 
in the Hyades 



HD no. 



V 



R-I 



mi 



[Fe/H] mi 



26756 


8.46 


0.35 


0.25 


0.01 


26767 


8.04 


0.31 


0.21 


0.36 


27771 


9.09 


0.39 


0.38 


0.26 


28099 


8.10 


0.32 


0.23 


0.01 


28258 


9.02 


0.43 


0.36 


-0.37 


28805 


8.66 


0.35 


0.29 


0.21 


28878 


9.38 


0.41 


0.42 


0.20 


28977 


9.65 


0.44 


0.45 


0.00 


29159 


9.37 


0.41 


0.40 


0.09 


30246 


8.31 


0.33 


0.24 


0.18 


30505 


8.98 


0.38 


0.37 


0.31 


32347 


8.98 


0.36 


0.31 


0.20 


284253 


9.14 


0.38 


0.35 


0.20 


284787 


9.05 


0.40 


0.37 


0.05 


285252 


9.00 


0.41 


0.43 


0.28 


285690 


9.56 


0.44 


0.52 


0.43 


285742 


10.26 


0.49 


0.59 


0.12 


285773 


8.94 


0.41 


0.36 


-0.12 


285830 


9.47 


0.44 


0.45 


0.06 


286789 


10.44 


0.52 


0.70 


0.47 


286929 


10.01 


0.51 


0.63 


0.15 



a; 
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Figure 4. The metallicities [Fe/H] mi for Hyades G and K dwarfs 
as a function of R — I. The dotted line represents the mean metal- 
licity value from the literature [Fe/H] = 0.11 (Taylor, 1994). Our 
mean value [Fe/H] = 0.16 is shown as a solid line. 



colour from the literature (Hauck and Mermilliod, 1998). A 
list of the Hyades G and K dwarfs, their magnitudes, B — V, 
R — I and mi colours and the metallicities calculated using 
eqn. (2) are shown in Table [l]. Taylor (1994) estimated the 
mean metallicity for the Hyades to be [Fe/H] = 0.11 ± 0.01. 
Our mean value is [Fe/H] = 0.16 ± 0.03. Fig || shows the 
relation between the derived metallicity and the R—I colour, 
and is found to be independent of colour, indicating that 
there are no residual temperature effects in the metallicity 
indicator. 



4.3 Comparison of the Stromgren and Geneva 
photometric metallicity calibrations 

In Fig ^ the relation between the metallicities calculated 
using the Geneva colour b\ and the Stromgren mi colour 
is shown. The separate metallicities are in good agreement 
with each other: the scatter around the fitted line is only 
0.15 dex, significantly less than the error of the individual 
metallicity estimates of 0.2 dex. This is probably because 
the b\ and mi colours are correlated, since they both mea- 
sure similar regions in the blue at approximately 4000 A. 
The metallicities are found to be in excellent agreement for 
[Fe/H] > —0.7, where most of the stars lie, differing by less 
than 0.01 dex. For the stars with metallicities [Fe/H] < —0.7, 
the scatter appears to increase and there may be some sys- 
tematic shift between the systems. In terms of the K dwarf 
problem, a possible systematic error at such low metallic- 
ity is not very significant. More low metallicity stars with 
good spectroscopic metallicity determinations would be of 
interest for testing the abundance indicators further. 




-2 -1.5 -1 -.5 .5 1 



[Fe/H] m , 

Figure 5. The relation between the metallicities based on the 
Stromgren m\ and Geneva b\ colours. The solid line shows a 
least square fit to the data (where the errors in both colours are 
taken into account). The dotted line shows the 1:1 relation. 

4.4 Metallicities from broadband photometry 

We have also developed a new metallicity indicator for K 
dwarfs using broadband BVI photometry. For this purpose 
we only choose calibration K dwarfs in the absolute magni- 
tude range 5.5 < My < 7.3. The calibration stars were se- 
lected from the Flynn and Morell (1997) sample, with B — V 
and V — I data being collated from Bessell (1990). The cal- 
ibrating sample is shown in Table El We recovered a simple 
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Table 2. K dwarfs used to calibrate a photmetric metallicity 
indicator based on B — V and V — I colours. Columns are the HD 
number, colours, absolute magnitude My, the spectroscopically 
measured metallicity [Fe/H] spe c and the photometric metallicity 
[Fe/H] BV i (Eqn. 3) 



HD no. 


B — V 


V — I 


My 


[Fc/H] sp ec 


[Fc/H]bvi 


4628 


0.890 


0.949 


6.376 


—0.40 


—0.29 


10700 


0.727 


0.802 


5.680 


—0.50 


—0.54 


13445 


0.812 


0.891 


5.930 


—0.28 


—0.50 


25329 


0.863 


1.093 


7.178 


— 1.70 


— 1.63 


26965 


0.820 


0.888 


5.916 


-0.30 


-0.41 


64090 


0.621 


0.865 


6.008 


-1.90 


-1.93 


72673 


0.780 


0.853 


5.953 


-0.35 


-0.48 


100623 


0.811 


0.870 


6.063 


-0.25 


-0.35 


103095 


0.754 


0.934 


6.611 


-1.40 


-1.33 


134439 


0.770 


0.959 


6.736 


-1.57 


-1.39 


134440 


0.850 


1.037 


7.077 


-1.52 


-1.31 


149661 


0.827 


0.832 


5.819 


0.00 


0.08 


192310 


0.878 


0.888 


6.002 


-0.05 


0.08 


209100 


1.056 


1.118 


6.893 


-0.10 


-0.17 


216803 


1.094 


1.169 


7.065 


-0.20 


-0.24 



relation between B 
follows: 



V, V — I colour and metallicity as 



[Fe/H] £ 



8.54 x (B - V) - 7.73 x (V - I) - 0.55 



(3) 



where we denote the metallicity index based on these 
colours by [Fe/HJevi- A comparison of the spectroscopic 
metallicities and this broadband photometric metallicity in- 
dicator is shown in Fig []. The scatter in the fit is ~ 0.1 dex. 
Note that this simple relation should not be used for stars 
with absolute magnitude My < 5.5, because the effects of 
stellar evolution begin to affect the colours in addition to 
the metallicity. 

The metallicity index is quite accurate, and is based 
on broadband colours only. As a consequence it is probably 
more useful than the indices based on the narrower band 
Stromgren and Geneva colours developed in the previous 
sections. A simple test of the calibration could be obtained 
by selecting low reddening open and globular clusters for 
which sufficiently deep main sequence photometry in B,V 
and I is available, since metallicities for these are known 
from spectroscopic studies. 

By applying this type of metallicity indicator to large 
photometric surveys (such as Sloan) it would be in principle 
possible to estimate the metallicity distribution of K dwarfs 
at different Galactocentric distances in the Milky Way's disk 
and even in the bulge region. This would represent an power- 
ful constraint for chemical evolution models. Moreover this 
metallicity indicator could also be used to study local group 
galaxies by utilizing deep two color data obtained with Hub- 
ble Space Telescope. 



4.5 Metallicities using K dwarf luminosity and 
colour 

We show in a companion paper (Kotoneva, Flynn and 
Jimenez, 2002) that stellar luminosity on the main sequence 
correlates very well with metallicity at a given colour. This 
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Figure 6. Relation between the metallicities using the broadband 
colours B — V and V — I for K dwarfs in the absolute magnitude 
range 5.5 < My < 7.3. The scatter around the 1:1 line is fs 0.1 
dex. 



was shown by measuring the displacement of stars AMy of 
K dwarfs from a fiducial isochrone in the My versus B — V 
plane (i.e. relative to the fiducial line at that colour). The 
fiducial isochrone comes from Jimenez, Flynn and Kotoneva 
(1998), has an age of 11 Gyr, solar metallicity, and was found 
empirically to be a good fit to solar metallicity K dwarfs. The 
metallicity determined by this method, [Fe/H]KF, is given 
by: 



[Fe/H] KF = 1.185 x AMy + 0.054. 



(4) 



This new metallicity indicator was developed on the ba- 
sis of the photometric derived metallicities described above. 
When we checked that [Fe/H]KF is consistent with the sam- 
ple of K dwarfs with spectroscopically determined metal- 
licities (and in which it is ultimately based, since the pho- 
tometric metallicities have been calibrated from these same 
dwarfs) the indicator turned out to be considerably more ac- 
curate than anticipated. The scatter between the [Fe/H]KF 
and the spectroscopically measured metallicities was « 0.08 
dex, i.e. little more than the intrinsic error in the spectro- 
scopic metallicities (0.05 dex). Although the spectroscopic 
sample is rather small, we consider that metallicities de- 
rived using this technique are greatly superior to photomet- 
rically derived metallicities. Note that the technique relies 
on accurate parallaxes (absolute magnitudes) being avail- 
able, and is at present applicable only to nearby stars (or 
stars for which an independent distance indicator is avail- 
able). Alternatively, the photometric techniques developed 
in this paper, although of lower precision, can be used on 
stars for which colours only are available. 

For all the stars in the sample we show the metallic- 
ity [Fe/H]KF in column 12 of Table [| The values are based 
on the Hipparcos B — V and My. The error in [Fc/H] is 
dominated by the colour error, which for the sample stars is 
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typically 0.025 mag, leading to a typical error in the metal- 
licities of 0.1 dex. 



5 K DWARFS IN THE HIPPARCOS CATALOG 
5.1 The data 

Part of the full dataset of 431 stars is shown in Table || 
The table shows the Hipparcos Input catalog (HIP) and HD 
numbers, visual magnitude V and absolute magnitude My, 
as computed from the Hipparcos parallax. The next five 
columns are the B — V colour, the mean values of R — I 
and mi and ci Stromgren colours and 6i from the Geneva 
catalog. The Stromgren based metallicity, [Fe/H] , Geneva 
based metallicity [Fe/H] 6 and the luminosity based metal- 
licity [Fe/H] KF follow. The three last columns show the U, V 
and W velocities in kms . The full dataset is available at 
the Strasbourg Data Center or from the authors. 



5.2 Raw metallicity distribution function for K 
dwarfs 

The normalised metallicity distribution function (MDF) for 
our basic sample of 431 K dwarfs is shown in the upper panel 
of Fig [j], and is based on the [Fe/H]KF metallicities. There 
are two studies of the metallicity distribution of K dwarfs 
in the literature with which we can compare our results. We 
compare the raw MDF, because this most closely matches 
the procedures which have been used to select K dwarf sam- 
ples in the past (rather than comparing to the more carefully 
selected sample of K dwarfs described in the next section) . 

Fig ^ also shows the metallicity distribution functions 
obtained by Marsakov and Shevelev (1988) and Rocha-Pinto 
and Maciel (1998b, RPM98). The Marsakov and Shevelev 
sample is an inhomogeneous compilation of metallicities 
from the literature. The Rocha-Pinto and Maciel stars were 
selected from the Third Gliese catalogue (Gliese et al, 1991), 
for which photometric data could be found from the litera- 
ture. The present sample is a near complete sample of the 
solar neighbourhood based on Hipparcos parallaxes and new 
photometric observations. The samples have rather different 
pedigrees, but they are all similar to the metallicity distri- 
bution for G dwarfs, being peaked at ~ —0.2 dex. 

The clear result in Fig ^ is that the metallicity distribu- 
tion we obtain here is broader than the other samples. We 
believe this can be attributed to the selection of our sample 
by absolute magnitude rather than by spectral type. We will 
argue this view in section 5.£ , and use what we regard as 
a superior procedure of selecting the stars by mass, rather 
than spectral type, colour or luminosity, as has been done 
in the past. 

Haywood (2001) has shown that the metallicity distri- 
bution of nearby G dwarfs for all samples in the literature 
peak at [Fe/H] « —0.2, as does the present K-dwarf sample 
and the RPM98 sample. Haywood (2001) however argues 
that the peak should be at [Fe/H] = 0.0, based on samples 
selected by colour and not spectral type, as has been done 
in the past. Our sample is selected by absolute magnitude 
rather than colour, since this more reliably selects the stars 
by mass, yet we still find the peak at its traditional location. 
A comparison of 101 stars in common between Haywood's 
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Figure 7. Normalised metallicity distribution functions (MDF) 
for K dwarf sampes in the literature and for this study. Lower 
panel : Rocha-Pinto and Maciel (1998b). Middle panel: Marsakov 
and Shevelev (1988). Upper panel: the 431 stars in our basic sam- 
ple, before any corrections. The three samples all show a similar 
spread in abundance and a peak at [Fe/H] ss —0.2. The present 
sample (upper panel) clearly shows a longer tail of super-solar 
metallicity stars; this is a conseque nce of the selection of the stars 
by absolute magnitude. In section 5.6 we develop a selection by 



stellar mass which substantially reduces the number of stars at 
super-solar metallicities. 



and our sample does, indeed, show an offset between our 
stars and his of 0.2 dex. Furthermore, a comparison with 
40 stars in common with Rocha-Pinto and Maciel (1998b) 
shows an offset of 0.15 dex (in the same sense as the Hay- 
wood comparison) . These two comparisons are not indepen- 
dent because the metallicities derived for the stars are par- 
tially based on Stromgren photometry in both cases. Closer 
comparison with the Haywood sample shows that the offset 
is a linear function of the absolute magnitude of the stars, 
rising from » 0.0 at M v = 5.5 to w 0.4 at M v = 7.3. Fig | 
shows the Haywood stars (the "long-lived dwarfs") in B — V 
colour versus [Fe/H]. The super-solar metallicity stars dom- 
inate preferentially among the cooler dwarfs, which suggests 
that there may be a systematic error in the metallicities as a 
function of colour (i.e. effective temperature) or some selec- 
tion effect. Subsequent to noting this trend, we found that it 
has already been commented upon by Reid (2002) . We note 
that restricting the Haywood sample to B — V < 0.7 would 
put the peak in the metallicity distribution at its traditional 
location (as we find here) at [Fe/H] ~ —0.2. We will later 
select our K dwarf sample by mass (section 5.6). The metal- 
licity distribution is not significantly altered if we divide the 
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Table 3. Part of the dataset showing identifications (HIP and HD numbers), visual V magnitude, absolute magnitude My, broadband 
colours B — V and R — I, Stromgren colours mi and ci, the Geneva colour 61, abundance estimates from various methods (see section 
5) and the space velocities U, V and W. The full table of 431 stars is available at the Strasbourg Data Center or from the authors. 
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Figure 8. Colour versus metallicity for the Haywood (2001) sam- 
ple of long-lived dwarfs. There is a systematic change on the 
metallicity distribution with stellar colour, with the super-solar 
metallicity stars dominating amongst the cooler spectral types. 

K dwarfs into subsamples by mass (see Fig |l3| ) within our 
adopted mass limits. 

5.3 K dwarf kinematics 

The basic sample consists of 431 K dwarf stars with metal- 
licity estimates accurate to » 0.1 dex. For 212 of the stars, 
radial velocities were found in the literature and space ve- 
locities U, V and W computed. Unfortunately, velocities 
are not available for all the stars; thus there is likely to be 
a small bias toward higher velocity stars in the literature 
sources, since high proper motion and metal weak stars are 
preferentially included in radial velocity programs. However, 
since about half the sample does have radial velocities, and 
the sample is furthermore dominated by thick disc and disc 
stars which have much smaller proper motions than halo 



stars, the bias in the measured velocity dispersions is likely 
to be quite small. A useful extension to this work would be 
to obtain velocities for all the stars. 

The velocities are shown as a function of metallicity in 
Fig |s| The velocity dispersions ajj , cry- and aw are shown 
as a function of metallicity in Fig |l^. The figures show 
the expected features of the solar neighbourhood: the thin 
disc with vertical velocity dispersion of about 20 km s for 
metallicities above approximately [Fe/H] = —0.5, the kine- 
matically hotter thick disc in the range —1 [Fe/H] ^ —0.5. 
There are very few stars with halo-like kinematics (high ve- 
locity dispersion and high asymmetric drift). This is as one 
would expect: for a local sample consisting initially of 431 K 
dwarfs, one would only expect of order 1 halo star, since the 
local disc:halo normalization is approximately 500:1 (Gould, 
Flynn and Bahcall, 1998). The apparent low metallicity stars 
are more likely the result of scatter from stars in the thick 
disc, since there is an observational scatter in the metallici- 
ties of 0.1 dex. 

5.4 Scale height /velocity correction 

Models of the chemical evolution of the local disc predict 
the metallicity distribution in a column through the disc, 
whereas the sample K dwarfs are drawn from a roughly 
spherical region centered on the Sun. The local sample is 
therefore biased toward stars of lower velocity dispersion 
since they will spend more time close to the Galactic mid- 
plane than older, faster moving stars. 

We correct for this by computing the velocity dispersion 
of the stars as a function of metallicity, and computing from 
this their vertical scale height using a realistic mass model 
of the Galactic disc (following Sommer-Larsen, 1991). 

The vertical velocity dispersions are tabulated in Table 
and shown in the upper panel of Fig [l^. In order to smooth 
out noise due to the small sample size, we have fit the ver- 
tical velocity dispersion linearly as a function of metallicity 
(shown as a solid line in the upper panel of Fig ^| We use 
the fit values of the velocity dispersion in what follows to 
determine the correction of the MDF for the scale height of 
the stars. 

The volume density of matter in main sequence dwarfs 
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Figure 9. Space motions U, V and W of the K-dwarfs as a 
function of the metallicity [Fe/H]KF- Note that the velocity scale 
is different in each panel. 
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Figure 10. Velocity dispersions crjy, cry and aw of the K dwarfs 
versus the metallicity [Fc/HJkf- The vertical velocity dispersion 
(jyy has been fit linearly with metallicity (upper panel) and is 
used to correct the raw metallicity distribution function for the 
effect of the scale heights of the stars. 



Table 4. Scale height correction for the metallicity distribution. 
The first column shows the metallicity ranges of the bins, the 
second column shows the observed velocity dispersion of the stars 
in the bin, ay^(obs), the third column the fitted vertical velocity 
dispersion, <rw (fit) (i.e. to the linear relation shown in the upper 
panel of Fighoh. The final column shows the ratio / of the column 
density to the local density of the K dwarfs computed in a model 
of the Galactic disc (normalised so that / = 1.0 at [Fe/H] = 
—0.02). The total mass density of the K dwarfs in the model is 
0.0043 M pc" 3 and the total column density (|z| < 2000.0 pc) 
is 3.5 M pc" 2 . 
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in the absolute magnitude range 5 < My < 8 is 0.0074 
M pc" 3 (Holmberg and Flynn, 2000, Table 1). The part 
of this which is represented by our sample K dwarfs (5.5 < 
M v < 7.3) is 0.0043 M pc" 3 (Holmberg, 2001, private 
communication). We have determined the vertical velocity 
dispersion for the K dwarfs as a function of metallicity, by 
sorting the stars by metallicity and dividing them into 11 
equal star number bins (so that each bin represents a lo- 
cal volume density of 0.0043/11 = 0.00039 M pc" 3 ). From 
this local density and the velocity dispersion of each bin, we 
then compute the total column density represented by each 
bin by integrating self-consistently the Poisson-Boltzmann 
equations in a model of the local Galactic disc (Holmberg 
and Flynn, 2000). A correction factor /, which is the ra- 
tio of the column density to the local density for each bin, 
normalised so that / = 1 in the bin closest to the solar abun- 
dance (i.e. at [Fe/H] = —0.02), is computed and is shown in 
Table §. 

5.5 Chromospheric activity corrections 

In this section we investigate the effects on the metallicity 
distribution of G or K dwarfs of chromospheric activity. 

Due to the Wilson-Bappu effect (Wilson and Bappu, 
1957; Wilson, 1976), photometrically derived metallicities 
differ somewhat from spectroscopic ones (RPM98, Rocha- 
Pinto and Maciel 1998a and references therein). The Wilson- 
Bappu effect in active chromospheres causes an emission line 
in the center of the stellar absorption lines, which leads to 
too low photometrically determined [Fe/H] metallicities. If 
a sample includes many active stars, the metallicity distri- 
bution would be biased towards metal poor stars. 

The effect on the present sample turns out to be quite 
small. Firstly, chromospheric activity is stronger among bi- 
nary stars. In the present sample the binaries have been 
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very effectively removed due to the high quality Hipparcos 
data for all the stars. The sample is therefore similar to 
RPM98, in which the binaries were also removed. Follow- 
ing RPM98, we can expect that about 30% of the stars the 
sample are chromosphericaly active after removing the bi- 
naries. The metallicity of these stars will be systematically 
incorrect. 

RPM98 have defined a correction to the metallicity dis- 
tribution, for metallicities determined by Stromgren pho- 
tometry, for G and K dwarfs. A similar analysis as in RPM98 
has been carried out for Geneva photometry (a medium 
band filter system quite similar to the Stromgren system), 
using those stars in our sample for which accurate spectro- 
scopic metallicities and measurements of the emission line 
width Rhk are available (Rocha-Pinto, private communica- 
tion). Although this yielded only 8 stars, the same trend was 
found for Geneva based metallicities as was found by RPM98 
for Stromgren based metallicities. Detailed computations 
show that the effect of this correction is to shift the mean 
of the metallicity histogram by ~ +0.05 dex, which can be 
easily understood since « 30% of the stars are thought to be 
chromospherically active and the mean correction to [Fe/H] 
for such stars is 0.15 dex (RPM98, their section 4). 

In the final sample of K dwarfs, we used metallicities 
based on the position of the stars in the Hipparcos colour 
magnitude (B — V versus My) diagram, rather than the 
purely colour based (intermediate band) metallicities above. 
Campbell (1984) has shown that the B — V colours of lower 
main sequence stars in the Hyades are affected by stellar 
activity, with typical changes in the B — V colour of order 
±0.02 mag; this would change the measured metallicities for 
the active stars by approximately ±0.1 dex. Unfortunately, 
without a detailed study, we cannot be sure if the effect 
would also yield a systematic error in the metallicities. 

In light of the above, we have decided not to make 
a chromospheric activity correction to the sample because 
we cannot quantify its effect well enough. The correction 
is likely to be small (w 0.05 dex), and probably shifts the 
position of the peak of the metallicity distribution rather 
than changing its shape. A major result of the paper, that 
the K dwarfs have a similar metallicity distribution to the 
G dwarfs, is therefore robust despite the uncertainties sur- 
rounding what chromospheric corrections should be applied, 
if any. 

5.6 K dwarf masses 

In the past, the metallicity distribution function (MDF) of 
the local disc has been determined for stars of a particu- 
lar spectral type or in some colour range. However, in the 
context of modelling the chemical evolution of the Galaxy 
and fitting it to the observed MDF, it is the range of stellar 
masses represented in the MDF which is of most interest. 
Selection of the stars within a particular range of masses 
permits a more direct comparison between the observations 
and theory. We describe in this section the use of isochrones 
to estimate masses for the K dwarfs. 

We show in Fig ^] the 5 Gyr Yonsei-Yale (hereafter Y 2 ) 
isochrones (Yi et al, 2001) for a range of metallicities. Also 
marked are the positions on each isochrone at which stars 
of masses 0.90, 0.85, . . . , 0.60 Mq lie. The two horizontal 
lines indicate the absolute magnitude cuts which were used 
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Figure 11. Yonsei-Yale (Y 2 ) isochrones for a range of metallic- 
ities. The positions of stars of various mass are marked with dif- 
ferent symbols. In the K dwarf region, 5.5 < My < 7.3 (dashed 
horizontal lines) there is a simple relation between mass, absolute 
magnitude My and colour B — V. 



to construct the basic sample, 5.5 < Mv < 7.3. Within 
the absolute magnitude range of interest, there is clearly 
a simple relation between mass and position in the colour- 
magnitude diagram. We have fit the mass M (in Mq) as a 
function of Mv and B — V as follows 



M = 0.322 x (B — V) — 0.169 x Mv + 1-551. 



(•>) 



The relation was found by fitting mass, colour and lumi- 
nosity in the Y 2 isochrones for ages ranging between 2 and 10 
Gyr. The internal precision in the relation is quite good, with 
a scatter of 0.03 Mq in the mass determinations. We have 
performed a similar analysis using the Padova isochrones 
for three metallicities and an assumed age of 5 Gyr, and 
find a very similar fitting relation for mass as a function of 
My and B — V as we did for the Y 2 isochrones. Masses ob- 
tained via the two isochrone sets follow a 1:1 relation closely, 
but there is an offset of 0.03 Mq , in the sense that the 
Padova masses are lower than the Y 2 masses. This difference 
is due to the differences in the adopted helium abundances. 
The Y 2 isochrones assume dY/dZ — 2 whereas the Padova 
isochrones assume dY/dZ — 2.25. The consequence of this is 
that the Padova isochrones have more Helium and so achieve 
the same luminosities as the Y 2 isochrones for slightly lower 
masses. Tests of the isochrones for Hipparcos K dwarfs in 
binaries and of known mass (Soederhjelm, 1999) would be 
of interest to better establish the scale and zero point of 
the mass calibration. For our purposes it is sufficient that 
the stars are marked by mass on a relative rather than an 
absolute scale, as is likely to be the case. 

Judging from Fig [111 our dwarfs typically have masses 
of 0.70-0.90 Mq. A more detailed plot of the K dwarf region 
is shown in Fig [l2|. It is clear that the absolute magnitude 
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Figure 12. Lines of constant mass for the Y 2 isochrones, plot- 
ted over the K dwarfs in our basic sample. The dwarfs were ini- 
tially selected by absolute magnitude 5.5 < My < 7.3 (shown 
by dashed horizontal lines). The absolute magnitude cuts bias 
the sample by including too many high metallicity stars at high 
mass and too many low metallicity stars at low mass. To avoid 
this we have restricted the sample to the mass range 0.75 < 
M/Mq < 0.83. 



cuts bias the sample (making it broader), by including too 
many metal rich stars at higher masses and too many metal 
poor stars at lower masses (see Fig ul upper panel). 

This is shown in detail in Fig O. We have marked on 
this plot the mass range in which the sample appears to 
be complete, i.e. between 0.75 and 0.83 Mq (on the scale 
calibrated to the Y 2 isochrones). Selecting stars within this 
interval reduces the sample from 431 to 220 K dwarfs. 

The metallicity distributions for the absolute magnitude 
limited sample (431 stars) and the mass restricted samples 
(220 stars) are shown in the lower panels of Fig [l4| The two 
distributions are quite similar; the main effect of the mass 
completeness restriction is to remove stars from the metal 
rich tail of the raw MDF. 

The upper panel of Fig ^ shows our final MDF for the 
220 mass selected K dwarfs after the velocity correct ion (i.e. 
by multiplying by the factor / described in section 5.4 and 
renormalising) . The velocity corrections slightly increase the 



relative fraction of metal weak stars and slightly decrease the 
relative fraction of metal rich stars, as expected. The MDF 
is shown in Table 



6 GALACTIC CHEMICAL EVOLUTION AND 
K DWARFS 

6.1 The theoretical model 

The model of Galactic chemical evolution we adopt here 
is that of Chiappini et al (1997) subsequently modified in 
Chiappini et al (2001), where a detailed description can be 
found. We review here the main ingredients of this model: 
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Figure 13. Metallicity [Fe/H] versus mass for our K dwarfs. The 
vertical lines show the limits within we judge the sample to be 
substantially complete by mass, i.e. between 0.75 and 0.83 Mq. 
The initial selection by absolute magnitude (to ensure that stellar 
evolutionary effects are negligible) causes the slanted edges to the 
left and right in the distribution of stars. 
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Figure 14. Normalised metallicity distribution functions (MDF) 
for our K dwarfs. The metallicities are based on [Fc/H]kf from 
Table ^. Lower panel : the open histogram is for the basic sam- 
ple of 431 K dwarfs, selected in the absolute magnitude window 
5.5 < My < 7.3. The filled histogram shows the raw normalised 
MDF for 220 K dwarfs selected to lie in the mass range 0.75 < 
M/Mq < 0.83 (c.f. Fig. [l]). The main effect of restricting the 
stars by mass is to remove some stars in the metal rich tail of 
the distribution. Upper panel : our final normalised metallicity 
distribution of K dwarfs in the vertical column at the Sun (i.e. 
including the velocity correction). This histogram is tabulated in 
Table |[ 
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Table 5. Tabulated metallicity distribution function for the final 
sample of 220 K dwarfs, corrected from local density to coloumn 
density (via the metallicity-velocity dispersion relation) and se- 
lected in the mass range 0.75 < M7M Q < 0.83. This is the MDF 
shown in the upperpanel of Fig [14. The metallicities come from 
[Fc/H]kf in Table fcil The error in the relative fractions is based 
on the Poisson sampling noise in the raw metallicity distribution. 
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• The model assumes that the halo + thick disc and the 
thin disc are formed during two different infall episodes. The 
thin disc does not form out of gas shed from the halo and 
the thick disc, but simply out of external gas. Such an inter- 
pretation is supported by recent dynamical and kinematical 
studies of stars in the outer Galactic halo by Sommer-Larsen 
et al (1997). Under these hypotheses, the infall rate is 



such as N and Fe, which are mostly produced by long- 
lived stars. 

• The prescription for the star formation rate (SFR) is: 

SFR oc 4>^ (7) 

where orot is the total surface mass density and a g is the 
surface gas density and ki = 1.5 and hi = 0.5. A threshold in 
the surface gas density is also assumed; when the gas density 
drops below this threshold the star formation stops. Two 
different values for the gas density threshold are assumed 
for the halo-thick disc (~ 4 Mqpc~ 2 ) and the thin-disc (~ 
7 M@pc~ 2 ) phases. The existence of such a threshold has 
been suggested by star formation studies (Kennicutt, 1989). 

• For the initial mass function (IMF) we adopt the pre- 
scriptions of Scalo (1986). 

• The contributions to the chemical enrichment from su- 
pernovae of different type (la, lb and II) as well as from stars 
dying as C-0 white dwarfs and contributing processed and 
unprocessed elements through stellar winds and the plane- 
tary nebula phase are taken into account in great detail (see 
Chiappini et al, 2001). 

• The adopted nucleosynthesis prescriptions are from: (a) 
van den Hoek and Groenewegen (1997) for low and interme- 
diate stellar masses and (b) Thielemann et al (1993) and 
Nomoto et al (1997) for SNe la. For the massive stars we 
considered the yields of Woosley and Weaver (1995) (case 
B) which include explosive nucleosynthesis. 



^) = x (X mf ) e-^ 

at (TTot(r,t G ) 

B(r) x 1 mS)l . , (6) 

where tt represents the time scale for the formation of 
the halo and the thick disc, and td{t) represents the time 
scale for disc formation, which is assumed to increase with 
Galactocentric distance. The best-fit model of Chiappini et 
al (2001) suggests that the timescale for the formation of 
the halo and thick disc is quite short and lie in the range 
~ 0.5 — 1.0 Gyr, whereas the timescale for the formation of 
the thin disc is quite long (~ 7.0 Gyr for the solar vicin- 
ity). This timescale for the thin disc ensures a very good 
fit of the new data on the G-dwarf metallicity distribution 
(Wyse and Gilmore, 1995; Rocha-Pinto and Maciel, 1996). 
A(r) and B(r) are derived by the condition of reproducing 
the present total surface mass density distribution in the so- 
lar vicinity. In other words, we integrate Eqn. ^ over time 
up to the present and normalise the left hand side by impos- 
ing that it is consistent with the disc's present total surface 
mass density. (Xi n f)i is the abundance of the element i in 
the infalling material, Tq the age of the Galaxy, assumed to 
be 14 Gyr and t ma x is the time of maximum gas accretion 
onto the disc coincident with the end of the halo-thick disc 
phase. 

• The Galactic thin disc is approximated by several in- 
dependent rings, 2 kpc wide, without exchange of matter 
between them. Continuous infall of gas ensures the tempo- 
ral increase of the surface mass density o Tot in each ring. 

• The instantaneous recycling approximation is relaxed. 
This is of fundamental importance in treating those isotopes, 



6.2 Observational and Cosmic scatter 

The model computes the metallicity distribution function 
(MDF) assuming no observational or cosmic scatter. The 
observational scatter in the metallicities for the K dwarfs 
is 0.1 dex (see section [l.E| ). This is quite small relative to 
the width of the observed MDF; the model can be easily 
convolved to take this into account when comparing to the 
data. 

Potentially more important is the cosmic, or intrinsic 
scatter, in the metallicities of co-eval stars. There is unfor- 
tunately no clear consensus presently regarding the amount 
of cosmic scatter. 

Edvardsson et al. (1993) found that the age- metallicity 
relation for F and G dwarf stars in the solar neighborhood 
shows a scatter of order 0.2 dex, which was larger than 
expected considering the uncertainties in metallicities and 
ages. This is of the order the width of the MDF, i.e. most 
of the MDF is accounted for by cosmic scatter, and the in- 
crease in the mean metallicity as a function of time plays a 
minor role. 

However, as discussed by Garnett and Kobulnicky 
(2000) the large scatter found for the stars in the solar 
vicinity is inconsistent with the abundance measurements 
in nearby spiral and irregular galaxies (e.g. Kobulnicky and 
Skillman, 1996) and in the local ISM (Meyer, Jura and 
Cardelli, 1998), which show that dispersions in ISM abun- 
dances are rather small on kiloparsec scales or less. By reana- 
lyzing the Edvardsson et al. sample together with Hipparcos 
parallaxes and new age estimates, Garnett and Kobulnicky 
(2000) found that the scatter in the age-metallicity relation 
depends on the distance to the stars in the sample. They 
concluded that the intrinsic dispersion in metallicity at fixed 



© 0000 RAS, MNRAS 000, 000-000 



K dwarfs and the chemical evolution of the Solar cylinder 13 



age is less than 0.15 dex for field stars in the solar neigh- 
bourhood, which is closer to the estimate of less than 0.1 
dex for Galactic open star clusters and the ISM (Twarog et 
al. 1997). 

Feltzing, Holmberg and Hurley (2001) argue also this 
view from a sample 5828 Hipparcos stars for which they de- 
rive metallicities and ages; they find a large scatter in metal- 
licity at any given age and even question the existance of the 
Age-Metallicity Relation itself. Rocha-Pinto et al (2000) ar- 
gue for a very small cosmic scatter, « 0.1 dex, based on the 
age-metallicity relation they present for nearby stars. 

As is well known (Pagel and Tautvaisiene, 1995) the 
metallicity distribution function is a very poor test of the 
Age-Metallicity relation, so we cannot resolve this issue here. 
We note that the AMR is less robust than the MDF. In 
the case of the MDF both quantities are directly observable 
(number of stars and metallicity), while in the case of the 
AMR the age scatter contributes significantly to the appar- 
ent metallicity scatter and is not easy to quantify. We will 
take the view that the amount of cosmic scatter is uncertain 
and we will keep it as a free parameter in the model fitting. 



6.3 Model comparison with the K dwarfs 

The metallicity distribution function for the model is com- 
pared to the data in Fig [l^. The model is a good fit. The 
data are shown in all panels by circles and with error bars 
based on the Poissonian sampling error (for the 220 star 
sample) . The solid line shows the model results based on the 
Woosley and Weaver (1995) and van den Hoek and Groe- 
newegen (1997) yields. This model is normalized to the value 
of the metallicity after 14.0 — 4.5 Gyr (i.e. when the Sun was 
born) . In the lower left panel of Fig |l5j the curves show the 
unconvolved model. In the remaining panels the model is 
convolved by a Gaussian representing the observational and 
cosmic scatter. The effect of the observational scatter alone 
is shown in the lower right panel (0.1 dex). The upper left 
panel shows the convolved model for a total scatter of 0.15 
dex, equivalent to 0.1 dex observational scatter and 0.11 dex 
of cosmic scatter, and the upper right panel shows a total 
scatter of 0.2 dex, equivalent to 0.1 dex observational scatter 
and 0.17 dex of cosmic scatter. Clearly, for this model, we 
cannot constrain the cosmic scatter beyond noting that it 
is less than 0.2 dex, and probably less than 0.15 dex. This 
is consistent with existing direct observational constraints 
(section 5.2 . We conclude that the agreement of the model 



prediction with the observed K-dwarf distribution is excel- 
lent, for any reasonable adopted cosmic scatter. An inde- 
pendent analysis of the width of the main sequence in the 
Hipparcos colour magnitude diagram by Girardi (2002, pri- 
vate communication) has also obtained the result that the 
cosmic abundance scatter is probably less than 0.2 dex. 



7 SUMMARY 

We have calibrated several photometric abundance indices 
for K dwarfs based on a sample of 34 G and K dwarfs with 
accurate, spectroscopically determined metallicities. Two of 
the indices use Cousins R — I photometry to estimate stellar 
effective temperature and the Geneva bi or Stromgren mi 
and ci colours. These indices give metallicity estimates of 
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Figure 15. Comparison of the model predictions and the metal- 
licity distribution for K dwarfs. In all panels the data are shown 
by circles. The lower left panel shows the two models, uncon- 
volved, as discussed in the text. The remaining panels show the 
model convolved with a Gaussian representing the amount of ob- 
servatinal scatter and intrinsic scatter. 



~ 0.2 dex accuracy. A third metallicity index uses B — V 
and V — I broadband colours. A fourth metallicity index, 
described in detail in a companion paper (Kotoneva, Flynn 
and Jimenez 2002, paper II), is based on absolute magnitude 
My and B — V colour. These latter two indices appear to 
provide metallicities of w 0.1 dex accuracy. 

For a set of newly acquired observations of K dwarfs, 
we use one of these indices to obtain the metallicity distri- 
bution of a near complete sample of K dwarfs in the So- 
lar neighbourhood drawn from the Hipparcos catalog. Care 
has been taken to remove the multiple stars from the sam- 
ple, for which metallicities cannot be measured accurately. 
Through isochrones we assign masses to the K dwarfs, and 
select those K dwarfs which fall into a mass window 0.75 < 
M/M0 < 0.83, within which the sample is near to complete. 
This yields a sample of 220 K dwarfs. 

The metallicity distribution of the 220 K dwarfs is 
strongly peaked near the solar metallicity, confirming the 
existence of the "G dwarf problem" amongst K dwarfs, 
as seen in several earlier studies (Marsakov and Shevelev 
(1988), Flynn and Morell (1997) and Rocha-Pinto and Ma- 
ciel (1998). We compare the metallicity distribution with 
Galactic chemical evolution models of Chiappini et al (2001). 
In the context of this model of the metallicity evolution, we 
find that the amount of cosmic scatter in the metallicities is 
small, not more than 0.15 dex. The model match the data 
well, indicating that the disc was formed via infall processes 
over an extended time scale of order 7 Gyr. 
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